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Abstract. A close relationship between the classical Hamilton- 
Jacobi theory and the kinematic reduction of control systems by 
decoupling vector fields is shown in this paper. The geometric in- 
terpretation of this relationship relies on new mathematical tech- 
niques for mechanics defined on a skew-symmetric algebroid. This 
geometric structure allows us to describe in a simplified way the 
mechanics of nonholonomic systems with both control and external 
forces. 
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1. Introduction 

The reduction of mechanical control systems to kinematic systems 
is very interesting and useful for solving control problems such as op- 
timal control problems |3] and for designing suitable control laws (see 
pi] . Chapter 8 in [1] and references therein). For instance, the plan- 
ning motion for the associated kinematic system determines trajecto- 
ries of the mechanical control system. Thus the methodologies to find 
these trajectories have been simplified because the kinematic reduction 
gives rise to a first-order control-linear system defined on the configu- 
ration manifold. Hence in general it is easier to solve or to analyze the 
kinematic system. If the mechanical control system is reducible to a 
kinematic one, then the controlled trajectories of this kinematic system 
under reparametrization define solutions of the original second-order 
problem on the phase space. An interesting particular case is the one 
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defined by kinematic reductions of order 1. Tliis kind of reductions de- 
fine a decoupling vector field. Unfortunately, tliere is not a systematic 
procedure for finding sucli kinematic reductions. 

Tlie pliilosopliy of kinematic reductions of order 1 seems, in a first 
approacli, quite similar to the standard Hamilton- Jacobi theory. This 
theory, that appeared with the dawn of analytical mechanics, is a valu- 
able tool for the exact integration of Hamilton's equations, for instance 
using the technique of separation of variables (see [1] and references 
therein). In many cases, the Hamilton- Jacobi theory allows us to sim- 
plify the integration of Hamilton's equations or, at least, to find some 
particular solutions. To be more precise, consider a configuration mani- 
fold Q and a hamiltonian function if : T*Q — )■ M. The Hamilton- Jacobi 
equation can be written as 



for some function W : Q M.. If we find such a function W, then the 
integration of the associated Hamilton's equations (with initial condi- 
tions along dW{Q)) is reduced to knowing the integral curves of a vector 
field Xfj^ on Q. This vector field is given by Xf}^ = TTr^QoXnodW e 
X{Q), where tt*q ■ T*Q — )■ Q is the canonical projection and Xh is the 
hamiltonian vector field associated to H . Hence, from the integration of 
a vector field on the configuration space it is possible to recover some of 
the solutions of the original hamiltonian system. Recent developments 
in Hamilton- Jacobi theory are described in |S1 El HSl [IH UHl 122] • Of 
course, the possible similarities with the theory of kinematic reductions 
are now clearer. 

One of the main objectives in our paper is to carefully study the 
underlying geometry of the kinematic reduction theory by showing the 
close relation with the classical Hamilton- Jacobi theory. Moreover, ad- 
vantage of recent developments in Hamilton- Jacobi theory for nonholo- 
nomic systems on skew-symmetric algebroids |2[ [H] even with external 
forces will be really useful to obtain a full novel theory of kinematic 
reduction for this type of systems. It is important to highlight this 
is not an arbitrary generalization since the mechanics on algebroids 
[HI ini HH |T71 [20] is particularly relevant for the class of Lagrangian sys- 
tems invariant under the action of a Lie group of symmetries including 
as a particular case nonholonomic dynamics (see [7] for a survey on the 
subject; see also [151 EHl 123] ) • 

The main results of this paper can be summarized in the following 
points: 
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• A description of nonholonomic mechanics in terms of the Levi- 
Civita connection associated to a fibered riemannian metric de- 
fined on the vector subbundle determined by the nonholonomic 
constraints. 

• A deduction of the Hamilton- Jacobi equation for nonholonomic 
systems in terms of the induced Levi-Civita connection. 

• An affine connection approach in presence of control forces. 

• A description of Hamilton- Jacobi equation with controls. 

• Relationship between Hamilton- Jacobi equation and kinematic 
reductions by decoupling vector fields. 

It is interesting to observe that our approach allows us to extend the 
theory of kinematic reduction to controlled system with symmetries as 
for instance, nonholonomic Lagrange-Poincare equations, etc. 

In the sequel, all the manifolds are real, second countable and 
The maps are assumed to be also C°°. Sum over all repeated indices is 
understood. 

2. Skew-symmetric algebroids 

In this section we introduce the notion of a skew-symmetric algebroid 
on a vector bundle td : D ^ M. It is known that this geometric 
structure covers many interesting cases in mechanics, as for instance, 
nonholonomic mechanics (see [E]). Similarly to the intrinsic definition 
of the Euler-Lagrange equations for a Lagrangian function L : TM — 
M obtained by the canonical structures on it (standard Lie bracket, 
exterior differential...), it is possible to determine the motion equations 
for a Lagrangian L : D — )■ M using the differential geometric structures 
naturally induced by the skew-symmetric algebroid structure. We will 
show that this generalization is quite useful in applications and clarifies 
the dynamics of systems with nonholonomic constraints. Let us first 
introduce the notion of a skew-symmetric algebroid. 

Definition 2.1. A skew -symmetric algebroid structure on the 

vector bundle '■ D ^ M is a ^.-bilinear bracket |-,-]_d '■ ^{td) x 
^{td) r(r£)) on the space T^td) of sections ofrn and a vector bundle 
morphism : D TM , so-called anchor map, such that: 

('^) [■) 'l-D is skew- symmetric, that is, 

lX,YjD = -lY,X]D, for X,Y e Tiro). 

(a) If we also denote by po '■ ^{td) — ^ X(M) the morphism of 
C'^{M) -modules induced by the anchor map then 

{XJY\d = flX,Y\D+pD{X){f)Y, forX,Ye T{td) and f G C^{M). 
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// the bracket |-, -Jd satisfies the Jacobi identity, we have that the 
pair {{■, -in, Pd) is a Lie algebroid structure on the vector bundle 
td : D M. 

If (|-, -IdjPd) is a skew-symmetric algebroid structure on the vector 
bundle T£, : D ^ M, then an almost differential of sections of 
A^td* being td* : D* ^ M the vector bundle projection of the dual 
bundle D* is defined as follows 

k 

(d^a)(Xo, Xi, . . . , Xk) = 5^(-l)Vz5(X,)(«(Xo, ...,X^,..., X^)) 

i=0 

+ J2i-iy+^a{lX„X,lD,Xo,Xr,...,Xi,...,Xj,...,X,) 

i<j 

for a e r(AVD.) and Xo,Xi, . . . ,Xk e T{td). 

In general (d^)^ 7^ 0. Indeed, (|-, -JdiPd) is a Lie algebroid structure 
on the vector bundle r^, : D ^ M if and only if (d^)^ = (see [HI HSl 
[201125] for more details about the Lie algebroids). 

Suppose that (x*) are local coordinates on M and that {e^} is a local 
basis of the space of sections T^td), then 

d 

The functions C^^, (pd)a ^ C°^(M) are called the local structure 
functions of the skew-symmetric algebroid on : -D — )■ M. 

If {e^} is the dual basis of {e^}, then 

where F G C°°(M) and k = KbC^ E T{td'). 

A po -admissible curve is a curve 7 : / C M — )■ D such that 

W = pDilit)) . 

Given X G r{To), the integral curves of the section X are those 
curves a : / C M — )■ M such that satisfy 

a = Pd{X) o a. 

That is, they are the integral curves of the associated vector field 
Pd(X) G X{M). If cr is an integral curve of X, then X o a is a po- 
admissible curve. Locally, the integral curves are characterized as the 
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solutions of the following system of differential equations 

= {pd)\X\x), 

where X = X^ca- 

Consider now the vector space over M 

H''{d^) = {f e C°°(M) I d^f = 0}. 

If M is connected and D is a transitive skew-symmetric algebroid, that 
is, 

Pj^{D^) = T.^M for all x G M, 

then 

H\d^) ~ M. (2.1) 

It is important to stress that condition (2.1 ) holds if the skew-symmetric 
algebroid has a connected base space and is completely nonholonomic, 
that is, 

Lie(~)(pz.(D))=T.M 
for all X G M . See [H] for more details. 



3. Bundle metrics on skew-symmetric algebroids and 

Newtonian systems 

3.1. The Levi-Civita connection. Let : D m D ^ M.he a, non- 
degenerate bundle metric on a skew-symmetric algebroid (D, |-, -Jd, pd)- 
Given this bundle metric we can construct a unique torsion-less connec- 
tion on D which is metric with respect to S (see [9] and references 
therein, for the standard case of Lie algebroids). The following construc- 
tion mimics the classical construction of the Levi-Civita connection for 
a riemannian metric on a differentiable manifold. 

We will denote by bgu : D ^ D* the vector bundle isomorphism 
induced by and by i^go : D* ^ D the inverse morphism. The 
bundle metric can be locally written as = {9^)abg^ ® ■ 

The Levi-Civita connection V^^ : ^{td) x T{t£)) — > ^{td) asso- 
ciated to the bundle metric is defined by the formula: 

2g^(vfy,z) = pz>(x)(g^(y,z)) + pz.(r)(g^(x,z)) 
-p^(z)(g^(x,y)) + g^(x, iz,Y]n) 
+g^(F,lz,x]^)-g^(z, iY,X}n) 

for X,Y,Z eT{TD). 

Alternatively, is determined by 

|X, YjD = VfY - Vl°X (symmetry) 

PDiX)i9^iY,Z)) = g^(vf r,Z) + g^(y, Vf Z) (metridty) , 
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These two properties allow to determine Christoffel symbols of the 
connection that satisfy 

More details about how to compute the Christoffel symbols are given 



in Section 5.3 if a S "^-orthogonal basis of r(rB) is taken. 

Additionally, we have the notion of derivative along an admissible 
curve. If7:JCM— )-Z}isa p^i-admissible curve and 

r(7) = {X: / C M ^ D|X is e°° and X{t) G D^^^t)'^ t e 1} 

is the set of sections along 7, then the induced covariant derivative 
V-yi r(7) — )■ r(7) can be defined as the mapping from X e r(7) to 
V-yX G r(7) with local expression 

'rl yC 1 

+ r^ABy^x^ 



v^x 



ec, (3.1) 



dt 

if 7(t) = {x\t),y^{t)) and X = X^ca- 

3.2. Geodesies. Given the bundle metric S^, a p^i-admissible curve 
7:JC]R— )-L)onL'is said to be a geodesic if 

V^7 = 0. 

If the local expression of 7 is 

7(t) = (x^(t),/(t)), 
then 7 is a geodesic if and only if 

The geodesies are just the integral curves of a vector field on D, called 
the geodesic spray ^go, whose local expresion is 

Note that if a G -D, then there exists a unique geodesic 7a : (— e, e) C 
D such that 7^(0) = a. If = r^, o 7^ : (-e, e) C M M is the 
base curve of 7^ then, since 7^ is p/j-admissible, we have that 

The associated symmetric product is defined as follows: 
(X : Y)sD = Vf r + Vf X , X, F G T{td) . 

The symmetric product on Riemannian manifolds is a fundamental tool 
in controllability results, kinematic reduction of mechanical systems 
and in the characterization of geodesic invariance of distributions (see 
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jlj). These results can be extended to our setting of skew-symmetric 
algebroids. 

Lemma 3.1. IfX,Ye T{td), then {X : r)^^ = [X'', [^g^,F^]], where 
X^ is the natural vertical lift of the section X of To- 

Proof. This is proved locally, similarly to the proof of this result on 
Riemannian manifolds in Lemma B.3 in [1]. 

Locally, X = X^ca and Y = Y^cb, then 



= X' 



d 



d 



dy 



A' 



Let us compute. 



c 



d 



dy 



c 



d_ 



Then, 



X' 



d 



d 



Y'iT% + V%^)y^ 



' dx 

d 
dy^ 



X\p 



d 
dy^' 

AdY^ 

dx"^ 



d)a- 



dx'' 



dX^ 
dx' 



+ y^iPDrs'^-^+x^y^inB+^'iA, 



d 



dy 



On the other hand, 

(X:r)gB = vfr + vfx = x^(p 



dY 



c 



dX 



c 



AvB-rC 



dx 
dX 



dY 



c 



dx' 



c 



+ ^''(p^)b^ + X^'Y^T^ + Fg^) j ec. 

This proves the equality since the vertical lift of this section is equal to 
the local expression of [X"^, [^g^, Y^]] we just computed above. □ 
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Now, we can extend the characterization of geodesically invariant dis- 
tributions already known on Riemannian manifolds [4J to skew-symmetric 
algebroids. A subbundle 2) of a skew-symmetric algebroid (D, |-, -J^i, p^)) 
is geodesically invariant if for any geodesic 7 : / C M — )• D with initial 
condition 7(0) G 2)(cr(0)), then 7(t) G D((T(t)) for any t G /, where 
o" = r£) o 7. 

Theorem 3.2. Let (D, |-, ■]/), p/j) be a skew-symmetric algebroid and 
Sd be a bundle metric over D. Let D be a subbundle of D and r^) = 
r£i|2). The following statements are equivalent: 

(i) D is geodesically invariant. 

(a) The restriction of the geodesic spray ^gu to D is tangent to D. 
(wi) IfX,Ye r(rs,), then \x,Y)^d G T{t^). 

Proof, {i) -v^ (u) The integral curves of ^gu are the geodesies. This 
proves the equivalence. 

{ii) =^ {Hi) Assume that X,Y E r(rj)). Then the restrictions of 
X^, to D are tangent to D. By hypothesis, (^313)1^, is tangent 
to D. Thus, [^go,F"^]|^ is tangent to D. As XJ^ is also tangent to 
D, we have that [X"^, [^go, y*"]] |^ is tangent to I). By Lemma 



3.1 



([X^[egz,,n])|^ = ((X: Y)lS\^. Thus (X: y)gz, G T{t^) because 
X^ restricted to D are tangent to D if and only if X G r(rj,). 

(Hi) =^ (i) Let {Xi, . . . ,Xr} = {Xa}a=i,...,r be a local basis for r(rj)). 
It can be extended into a local basis for r{To), which is given by 

{Xi, . . . , Xr, Xr+l, . . . , Xn} = {Xq, X^}, 

for a = 1, . . . , r; a = r + 1, . . . ,n. Hence, 

Vi>, = r^,X, + r°X„, for a,be{l,...,r}. 

By assumption, (Xq: X;,)go G T{tx,), then F"^ -|- F^^ = for all a, 6 G 
{1, . . . , r}, a E {r + 1, . . . ,n} because 

{Xa : X;,)gD = V^^Xfe + V^^^X^ = (F^^ + Fg„)Xc + (F"f, + F^„)X„. 

Now, let 7: /—> D be a geodesic such that 7(0) G 2)(cr(0)) where 
a = td o •y. Suppose that 

r n n 

7(t) = 5^«"(t)X,(a(t)) + ^ «"(t)X„(a(t)) = ^«^(t)X,(a(t)), 

a=l a=r+l i=l 

then u"(0) = for all a; G {r -|- 1, . . . ,n}. As 7 is a geodesic, 

— + Ff,.(a(t))nXtK(t) X,(a(t)). 
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Then 

— + (rf^. o a)wV = 0, \fke{l,...,n}. (3.2) 

As + = for all a, 6 G {1, . . . , r}, a E {r + 1, . . . , n}, we have 
that M'^(t) = for allt G / is a solution of the differential equation (3.2) 
for a = r + 1, . . . , n: 

= (Coa)«V = 0. 

As the initial condition for is u°'{0) = for all a G {r + 1, . . . , n}, it 
follows that u^lt) = is the unique solution for the differential equa- 
tions in (3.2). Thus, 



7 



(t) = J2 u''{t)Xa{a{t)) G 2)(a(t)), V t G / C 



a=l 



Hence D is geodesically invariant. □ 



3.3. Newtonian systems. Given a bundle map "J : D ^ D (that is, 
To = To) we define a newtonian system as the triple {D, Q^, 5"). 
This newtonian system induces the system of differential equations: 

VS",)7(t) = Hlit)), tel, (3.3) 
where the solutions are curves 7 : J C M — )■ D which are p£)-admissible. 

Given local coordinates {x\y^) associated with the basis {e^} for 
sections of D, Equations (3.3) can be written as 

where 3^{x\y^) = {x\J'^{x^,y^)). 

Observe that Equations (3.4) are the equations of the integral curves 
of a vector field ^go^g- on D. Locally, this vector field is given by 

Remark 3.3. The map 5" could be given by a section F G ^{to) such 
that 3^ = F o To- An interesting particular case is when we have a 
potential function V : M ^M. and F is the section — gradgoV^ G T{to) 
given by 

g^(gradgi,\/,X) = Pd{X){V) = dV{po{X)), for every X G T{to). 
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In particular, the solutions of this newtonian system (D, S^, = 
—grad^DV) are equivalent to the solutions of the Euler-Lagrange equa- 
tions on a skew-symmetric algebroid with Lagrangian L : D — > M: 

L{v) = l9''{v,v)-V{To{v)). (3.5) 

Therefore, these solutions are p£)-admissible curves 7 : / — > D such 
that 

V^^l^lit) + gradgz,F(rz5(7(t))) = 0. (3.6) 
Locally, those solutions satisfy 

y - abv y -{b ) [pd)bq-[ ^ 

where (3^)^^ are the entries of the inverse matrix of {{9^)ab)- In this 
case, the vector field p on D is given by ^go y = ^go -|- {grad^oVy , 
where {grad^oVY is the vertical lift to D of the section grad^oV G 
r(r^). o 

Example 3.4. If D = TM, I,]/? = [, ] the standard Lie bracket on 
M, pu = Id-TM and is a riemannian metric on M, then Equations 



(3.6) are the classical Euler-Lagrange equations for the mechanical 
lagrangian L : TM — y M. 

Example 3.5. Given a regular distribution D on TM and a riemann- 
ian metric S"^*^, we consider the riemannian orthogonal decomposition 
TM = D Q) and the associated orthogonal projectors P : TM D 
and Q : TM D^, see |2l [H]. Denote also by : ^ ^ TM 
the canonical inclusion. We induce by restriction a bundle metric 
X ]\.f D — } ]R and an skew-symmetric algebroid structure on 
D as follows: 

where X, F G r(r£)). Note that in this example, X, Y are vector fields 
on M taking values on D. Moreover, the Levi-Civita connection 
coincides with the constrained connection V^F: = P(V^y) defined, 
for instance, in |1] if V'^ is restricted to r(r£)). 

For this particular skew-symmetric algebroid structure. Equations 



(3.6) correspond with the equations of the nonholonomic system de- 



termined by the constraints induced by the distribution D and the 



mechanical lagrangian (3.5). These equations are also called in the 
literature Lagrange-D'Alembert's equations. 

Consider a basis of S'^-orthogonal vector fields {X^jXa}, 1 < ^ < 
m = rankD, m + l<a<n = dimM, adapted to the decomposi- 
tion TM = D Q) D^. In other words, D^^ = span {Xa{x)} and = 
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span {Xa{x)}. Observe that for the induced coordinates {x\y^,y°') on 
TM the nonholonomic constraints are rewritten as = 0, m + 1 < 
a < n. That is, the induced coordinates on D are given by {x\y"^). 
Therefore, the skew-symmetric algebroid structure induced on the vec- 
tor subbundle D — > M is locally described by: 

{Xa, Xb]d = P[Xa, Xb] = P{CabXc + ^ab^p) = ^ab^c, 
Pd{Xa) — Xa , 

where Xa = ipDYA£i, l<A<m. 

The Lagrange-D'Alembert's equations are 

= {pDTAy^, 

f = -T^y^y'-i^TiPDYB^. (3.7) 

Example 3.6. Our theory is not only restricted to lagrangian systems 
defined on the tangent bundle TM or nonholonomic systems determined 
by a regular distribution on TM. The techniques described in this paper 
by means of skew-symmetric algebroids are general enough to cover the 
most important cases of reduction of mechanical systems subjected or 
not to nonholonomic constraints. 

As a particular example, we include in our analysis the case of Lie 
algebras of finite dimension (it is clear that is a Lie algebroid over 
a single point). Now, suppose that (/,0) is a nonholonomic Lagrangian 
system on g, where Z : g — >■ R is a Lagrangian function defined by 
/(^) = J : — )■ 0* is a symmetric positive definite inertia 

operator and i) is a vector subspace of 0. We have the orthogonal 
decomposition 

= oeo^, 

where 0"'" = G | = G O}. Take now an adapted basis to 

this decomposition {6^,6^} where t) = span {e^} and D"*" = span {cq}. 
Then, the Euler-Poincare-Suslov equations for (/,()) are 

y — ^ Asy y i 

where {|/^, y"} are the global coordinates on induced by the basis 

{eA,e„}. 

Example 3.7. Similarly, more involved situations can be recovered 
using our techniques. For instance, nonholonomic systems on Atiyah 
algebroids associated with principal G-bundles. For the sake of simphc- 
ity, we will consider the particular case when the principal G-bundle is 
trivial. In such a case, the Atiyah algebroid is a vector bundle of the 
form 

ta- A = QX TM — > M, 
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where g is the Lie algebra of the Lie group G and M is a smooth 
manifold. The Lie bracket of the space r(r^) is characterized by the 
following condition 

for ^, ^' G g and X,X' G X(M). The anchor map pA is the canonical 
projection onto the second factor. 

Suppose now that D is a vector subbundle of A over M of constant 
rank (the constraint bundle) such that 

M3x-^ Dv{x) : = D{x) H (g x {Ot^m}) Cqx T,M 

is a vector subbundle of A. Then we can choose a local basis {^a}i<a<r 
of r^Toy), with C,a '■ U M — )■ smooth maps, and a local basis 
{Xa} = iVa, Y^)} of r{TD), With : f/ C M and r„ G X{U). 

Moreover, if (x*) are local coordinates on f/ C M and = Y^^, 
the Lagrange-D'Alembert-Poincare equations are: 

= Y^^y'^, 

r = -T%y\ 

where {x^,y'^,y°') are the corresponding local coordinates on D. 

Note that in the particular case when M is a single point, we recover 
the Euler-Poincare Suslov equations. 

Geometric interpretations of nonholonomic LR systems or non- 
holonomic systems with semidirect product symmetry may also 
be given using skew-symmetric algebroids deduced from Lie algebroids 
(see |8] for more details). 



/ nD\ca 

- ^>^.(s ) , 



4. Hamilton- Jacob: equation 



The next result is a direct consequence of Equations (3.6). See [2] for 
an extension of this result (in a hamiltonian context) for many differ- 
ent types of mechanical systems (nonholonomic dynamics, dissipative 
systems...). 

Proposition 4.1. Let {D, | ,]d,Pd) be a skew- symmetric algebroid and 
consider a newtonian system determined by {D,S^,3^)- Take an arbi- 
trary section X G ^{td) then, the following conditions on X are equiv- 
alent: 
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(i) If a : I — y M is an integral curve of the vector field Pd{X) that 
is, 

a{t) = pn{X){a{t)), (4.1) 
then J = X o a : I — > D is a solution of 

(a) X satisfies 



Proof. Let a : / — )■ M be an integral curve of the vector field pd{X). 
Then, the result follows immediately from the fact that 

Vjl(,)(Xoa(t)) = (VrX)oa(t). 

□ 

This result is analogous to the Hamilton- Jacobi theory already de- 
scribed in a Lagrangian framework for a free mechanical system in [5] 
and for a nonholonomic one in [0], but adapted to skew-symmetric al- 
gebroid structures. 

Remark 4.2. An analogous result can be written replacing the section 
X by an entire distribution spanned by sections. o 

Remark 4.3. In standard riemannian geometry (that is, D = TM 
equipped with the standard Lie bracket and without external forces 
5" = 0) the vector fields X satisfying X = are called auto-parallel 
vector fields and are obviously connected with the solutions of the geo- 
desic equations. o 



From now on, we only consider mechanical problems given by (D, S 



D 



V) as described in Remark 3.3 Specializing Proposition 4.1 to this 
kind of problems and to vector fields X verifying an extra condition 
ij>cd^(bg(X)) = we obtain a new expression of this Proposition. In- 
deed, the following Theorem can be compared with the classical expres- 
sion of the Hamilton- Jacobi equation proposed in [ll] for the hamilton- 
ian function h : D* — > M: 

h{K) = g^*(/€, ft:) + r(r^. (/€)), 

where k e D* and S^* : D* 

X M D* — ^ is the induced bundle metric 

on the dual bundle. 

Theorem 4.4. Let (-D, [ ,}d,Pd) be a skew- symmetric algehroid and 
consider a mechanical problem determined by [D,'^^ ,V). Take a sec- 
tion X G '^{td) such that ijs:d'^(bgo(X)) = 0. Under this hypothesis, 
the following conditions are equivalent: 
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(i) If a : I — y M is an integral curve of the vector field po{X), 
that is, 

a{t) = pn{X){a{t)), 
then 7 = X o (T : / — > D is a solution of 

+ grad^uV{rD{l{t))) = 0. 

(ii) X satisfies the Hamilton- Jacobi differential equation 

d^Qg^(X,X) + y^ =0. (4.2) 

//, additionally, the skew- symmetric algebroid (D, |-, ]/), p/)) is com- 
pletely nonholonomic and M is connected or if H^{d^) ~ M, then Equa- 
tion (4.2) is equivalent to 

^g^(X,X) + 1/ = constant. 

Proof. First, we study how to express the condition ixd^{\>c,D{X)) = 
in terms of the Levi-Civita connection associated to S^. Let Y G T^td), 

= d^(bg.(x))(x,r) 
= pz5(x)(g^(x,F)) -pz5(r)(g^(x,x)) - g^(x, |x,ri^) 
= g^(vr X, Y) + g^(vr F, X) - g^(vrx, x) - g^(vr x, x) 

-g^(x,vrr-vrx) 
= g^(vfx,F)-g^(vrx,x). 

Therefore, the condition ixd^{\>5D{X)) = is alternatively written as 
g^(Vf X, Y) = g^(VS°X, X) for every Y e T{td). (4.3) 



We only need to check that both condition (ii) in Proposition 4.1 



and Theorem AA are equivalent. If we examine the Hamilton- Jacobi 
differential equation 

^ig^(x,x) + v^ =0, 

then for any Y G T{td) and X G T{td) satisfying ixd^(t'gi5(X)) = 
we have that 

= d^(^g^(x,x) + v)(F) 

= lpD{Y){9''{X,X))+pn{Y){V) 

= g^(vfx,x) + pz,(F)(v) 

= g^(Vf X, Y) + g^(gradgz,r, Y). 
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In the last equality we have used condition (4.3). Therefore Equation 



(4.2) is written as 

Vf'x + gradgi,l^ = 0. 
(See also [21 [H]). □ 



From Theorem 4.4[ it is clear that we need to find sections X satis- 
fying ixd'^(bgD(X)) = 0. The most simple candidate to be a solution 
is X = gradgo/ with / : M -> M. Note that if {d^)^{f) = 0, then 
ixd^{\>S'^iX)) = 0. After some straightforward calculations, this con- 
dition is equivalent to: 

[pn{Y),pn{Z)]{f) = {pnlY,Z}n){f) . 

This condition is always true if the bracket [tId satisfies the Jacobi 
identity, that is, if the pair (|-, -JdjPd) is a Lie algebroid structure 
on the vector bundle To '■ D ^ M (see 



5. Mechanical control systems and kinematic reductions 

Assume that the newtonian system determined by (Z), S"^, 3^) also 
contains some input forces. We model this set of input forces by a 
vector subbundle X)^"^^ of D* . Locally, 2)'^'^) = span {6*^, . . . ,6*^^}, where 
6*' G r(r£).), 1 < / < A;. Denote by !D(c) the vector subbundle of D 
defined by D(c) = tlgi5(!D*^^^). Therefore, locally I'(c) = span {Yi, . . . , Y^} 
where Yi = Hgo^', Yi G T{td), 1 < I < k. The vector fields Fi, . . . , 
are called the control sections or input sections. 

The equations of motion for a newtonian system with input sections 
are as follows 

V?",)7(t)-3^(7(t))e%)(7(t)), VtG/CM, (5.1) 
where 7 : / — > is a p^-admissible curve. 



In terms of the control sections. Equation (5.1) can be rewritten as 
follows: 

k 

V?w7(t) - S^ilit)) = Y,u\t)Yiirom)) (5.2) 

1=1 

for some u : / C M — > M*^, playing the role of controls. The correspond- 
ing local equations are 

Note that if F G V{tt)^^^) then the integral curves of the vector field 
^gD gr -|- Y^ on D are solutions of the previous equations, where Y^ is 
the vertical lift of the section Y . 



16 



BARBERO, DE LEON, MARRERO, MARTIN DE DIEGO, AND MUNOZ 



Remark 5.1. It is also possible to study the more realistic case when 
f/ is a proper subset of M*^. In this case, the controls u take value in 
a proper set of M'^ (i.e., not all the linear combinations of controls are 
allowed). Our procedure can be adapted to that particular control set. 
But for geometrical clarity in this paper we only consider the control 
distribution ©(c)- o 

Definition 5.2. The 4-tuple {D,S^ ,3^,'^{c)) is called a mechanical 
control system on a skew- symmetric algebroid. 

Consider the orthogonal decomposition D = !D(c) © 2)^^ induced by 
the bundle metric 9^ , with associated orthogonal projectors P(c) : D — )• 
and g(,) : D ^ T)f^y 

The following proposition is a direct consequence of the definition of 
mechanical control systems. 

Proposition 5.3. An admissible curve 'y : I ^ D is solution of Equa- 



tion (5.1) if and only if ^ : I ^ D satisfies 

Qic) (V?"i)7(t) - Hlit))) =0 V t G / C M. 

5.1. Kinematic reduction of mechanical control systems on a 
skew-symmetric algebroid. Now we introduce the notion of a kine- 
matic reduction (see [4J). 

Given an skew-symmetric algebroid structure on the vector 
bundle td : D M, we define a driftless system as the set (M, 2), U), 
where 2) is a vector subbundle of D locally spanned by {Xi, . . . , X^'}, 
with Xq, G r(r£i), 1 < a < k', and t/ C M'^ is the set of admissible 

controls. For a section X = Ylti=i^°'-^a £ ^{'Td), remember that an 
integral curve of X is a curve a: / C M — )■ M such that 

k' 

d(t) = J]n"(t)p^(X«)(a(t)), (5.3) 
«=i 

where iu\t), . . . ,#'(t)) = {u\a{t)), . . . ,u''' {a{t))) E U for all t. It can 
also be said that the pair (cr, u) is a solution to the driftless system. 

Observe that for each pair {a(t),u{t)) we have the curve 7 : / — 
D <ZD defined by 

k' 

l{t) = Y.u-{t)X^{a{t)). 
In the sequel, we will denote by the restriction {t£,) . 
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Definition 5.4 (Kinematic reduction). Let (D, 9^, 3^, 2)(c)) be a me- 
chanical control system. A driftless system {M,T),U) is called a kine- 
matic reduction o/ (D, S^, 3", D(c)) if for every solution {a{t),u(t)) 



T!:=^u-{t)X^{a{t)). 

The rank of a kinematic reduction is the rank of the distribution D. 
Rank-one kinematic reductions are particularly interesting. A section 
X of r(r£)) is called a decoupling section if the rank-one kinematic 
system induced by 2) = span {X} is a kinematic reduction. 

Definition 5.5 (Kinematic controllability). A mechanical control sys- 
tem (Z), S^, 5", !D(c)) is kinematically controllable if it possesses de- 
coupling sections {Xi, . . . ,Xk'} whose involutive closure has maximum 
rank. 

When a system is kinematically controllable, motion planning is pos- 
sible by using concatenations of integral curves of the decoupling vector 
fields. Those curves must be reparametrized in such a way that each 
segment begins and ends with zero velocity, see |1] for more details. 

We have the following adaptation of the results in Section 4. 

Proposition 5.6. Let {D, | ,}d,Pd) be a skew- symmetric algebroid and 
consider a mechanical control problem determined by (D, S^, 5", D(c)). 
Consider a driftless system (M, D, U). For all X G r(rj)), the following 
conditions are equivalent: 

(i) if a : I — > M is an integral curve of p£,{X), that is, 

a{t)=pD{X){a{t)), 

then 7 = X o (T : / — y D is an admissible curve solution of 

(tt) Q{c) (vf X-3^ox) =0. 

Under extra assumptions, we have an alternative way to write condi- 



of (5.3) there exists a pair {'y{t),u{t)) solution of (5.2), where 7(t) 



tion (ii) in Proposition 5.6 This provides us with a new characterization 
of kinematic reductions in terms of the affine connection of the given 
mechanical control system. (See also |H l2T]). 

Proposition 5.7. // (5(c) (?") = and D = span {Xi, . . . ,Xk>}, then 
the following conditions are equivalent: 

(%) For allX&T{T^), Q(e) (vf X-3^ox) = Q(,) (vf x) = 0. 
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(ii) For a//a,/3,7e {l,...,/c'}, = 0, Q(e)((X^ : X^)g^) = 

0. In other words, for all a, (3,'-]' e {1, . . . , k'}, G V{td^^^^) 
and {X^. X^)s^ G Tir^s^J. 
(Hi) {M,D,U) is a kinematic reduction o/ (D, S'^, 3^, 2){c))- 

Proof, (i) =^ (ii) Observe that for each a = 1, . . . ,k' we have that 
fX^ G r(rx,) for all / G C°°(M). If (i) holds, we have 

= Qic) {"^fxJXa) = f'Q(c) (V£X,) + fpniXM)Qic)iX^) 

for every / G C°^(M). Then, taking suitable functions /, we obtain 
Q{c)iXa) = and (5(c)(V^^Xq,) = 0. Now, using the polarization 
identity we have that 

{Xa : Xfj)^^ = - (Vx^^x^i^o, + X/s) - Vx^Xa, - Vx^Xf. 

Therefore, Q(c){{Xa ■ ^/3)sb) = 0- 

(ii) =^ (iii) By definition D (1 D. As (ii) is true, then statement (ii) 
in Proposition |5.6| is satisfied for every X G T{t£)). Hence statement 
(i) in Proposition 5.6 is true and we can conclude that (M, D , U) is a 
kinematic reduction of {D, S'^, 5", I'(c)) according to Definition 



5.4 



(iii) =^ (i) As (M, 2), U) is a kinematic reduction of {D, S , 3^), state- 
ment (i) in Proposition |5.6 is satisfied. Hence the result follows. □ 



A straightforward corollary of Proposition 5.7 is the following one: 

Corollary 5.8. If Q(c) {^) = and D = span {Xi}, then the following 
conditions are equivalent: 

(i) Qic) (v|>i - 3^ o Xi) = g(,) (vi'Xi) = 0. 

(a) (M, D, U) is a rank-one kinematic reduction of {D, S^, 3", ©(c))- 

In the sequel assume that 5" comes from a potential function V : M ^ 
M, that is, 3" = — gradguT^ o td. Then Theorem 4.4 can be adapted to 
control systems as follows. 

Theorem 5.9. Let (-D, | ,}d,Pd) be a skew- symmetric algebroid and 
consider a control problem determined by {D, Q^, V, D(c))- 

Take a section X G ^{td) such that 2xd'^(bgi?(X))(y) = for all Y G 
T{t^± ). Under this hypothesis, the following conditions are equivalent: 

(i) If a : I — M is an integral curve of pd{X), that is, 

a{t) = pn{X){a{t)), 
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then 7 = X o (T : / — y D is an admissible curve solution of 

Qic){V%^l{t) + gradsnViTDiiit)))) = 0- 
(a) X satisfies the Hamilton- Jacobi differential equation 

Qg^(X,X) + (Y) = for all Y G Tir^^^J . (5.4) 



Proof. Following similar arguments that in the proof of Theorem AA 
it is easy to prove that the condition ixd^(\>sD{X))(Y) = for all 
Y G r(r23± ) is equivalent to 

g^(vrx,y) = g^(vrx,x). 

Finally, observe now that 

= d^Qg^(x,x) + F^(r) 

= g^(vfx,x) + p^(r)(i^) 

= g^(vfx,y) + g^(gradgz,v,F), 

The last expression is equivalent to the equation 
g(,)(Vfx + gradgz,V)=0. 

□ 

Hence, we have just extended the notion of decoupling sections for 



mechanical control systems with nonzero potential since Theorem 5.9 
gives sufficient and necessary conditions to have a kinematic reduction 
of such a mechanical control system. 

Remark 5.10. Note that the condition zxd^(bgD(X))(F) = in the 
hypothesis of Theorem 5.9 is C°°(M)-linear in F G "^fc)- Hence only for 



a basis of vector fields in 2)^^ the condition must be checked. However 
that same condition is not C°°(M)-linear in X G r(r£)). o 

This theorem plays a key role to define in the future the geometric 
notion of motion planning for mechanical control systems with non-zero 
potentials, only known so far for zero potentials [3]. In that sense it will 
be useful to have some notion of reparametrization of integral curves of 
sections X G T{td). 

Proposition 5.11. Let X be a nonzero section in r(rDj^j) such that 
for all Y G T{t<j^± ), ixd^(\>QD{X)){Y) = 0. For a nonzero function f 

m e°°(M), the section fX G T{t'o^^^) satisfies ifxd^{^^o{fX)){Y) = 
for all Y G FfrDX ) if and only if PDiY)(f) = for all Y G Ffr^jx ). 
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Proof. Let us rewrite the condition ifxd^{\>QD{fX))(Y) = as follows 

= g^(vsj/x,y)-g^(vf /x,/x) 
- /pz,(r)(/)g^(x,x) = -/pz5(r)(/)g^(x,x) 

because X G r(rDj^j). From here, the equivalence is straightforward. 

□ 



As a consequence of Proposition 5.11 Theorem 5.9 can also be writ- 
ten for such a fX. 



From Theorem 5.9, we establish a connection between decoupling 
sections in the sense defined in [Ij and the solutions to Hamilton- Jacobi 
differential equation. 

Corollary 5.12. Let {D, | ,Jd,Pd) be a skew- symmetric algebroid and 
consider a control problem determined by {D, g^, V, 2)(c)) where V sat- 
isfies that gradgcV G T{txi^^^)- If a section X G r(r23^^j) is a decoupling 
section, then the following conditions are equivalent: 

(i) zxd^(bgD(X)))(r) = for all Y G T{t^i^J, 

fti) = g^(vf X, X) for all Y eT{T^ 



(Hi) X satisfies the Hamilton- Jacobi differential equation 

d^ Qg^(X,X) + (Y) = for all Y G Tir^.J . 

Proof The condition ixd^(bgD(X))(F) = for all Y G r{T^± ) can be 

g^(vfx,r) = g^(vfx,x) 



rewritten as 



for all y G r(r. 



As X is a decoupling section, both X and X G r(rDj^j) by Corol- 
lary 5.8 Hence the left-hand side of the above equality is zero and [(I) 
is equivalent to |(ii 



Condition 
0= d^ 



can be rewritten as follows: 



^g^(X,X) + (F) = lpoiY)9''{X,X) + pn{Y){V) 



= g^(vf x,x) + g^(gradgz,v^,y), 

for all Y G T{t<j)± ). Hence the equivalence between 
because of the hypothesis for V. 



and 



is clear 

□ 
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Note that the property of being a decouphng section is preserved by 
e°°(M)-multiphcation. However, the condition zxd^(bgD( X))(y ) = 

has al- 



for all Y G r(rD± ) is not C°°(M)-linear on X. Proposition 



5.11 



ready characterized those functions / such that i/xd (bgD(/X^ 
for all Y G T(t^± ). 



It might also be the case that a decoupling section X G r(r23 



does not satisfy |(ii) in Corollary 5.12, but there might exist functions 



/ G e°^(M) such that fX satisfies (ii) Let F G r(r^ 



g^(vr/x,/x) = f9''iPDiY)if)x + fvfx,x) 



■>D 



fpDiY)im^iX,X) + mVfX,X). (5.5) 



Condition (ii) in Corollary 5.12 would be satisfied if the following partial 
differential equation for / has solutions: 

fpniY){m''{X, X) + /2S^(Vf X, X) 







for all Y G Vir. 



The chances to find a solution depend on the 



particular examples under study, see Section |5.3[ 

5.2. Maximally reducible systems. If {M,D,U) is a kinematic re- 
duction of (D, S"^, 3^, D(c)), then any solution of (5.3) can be followed 



by a solution of (5.1). In this section we consider when the converse 



is also possible in such a way that we can talk about "equivalence" of 
controlled trajectories as mentioned in the introduction. 

Definition 5.13 (Maximal reducibility). A mechanical control system 
{D , , 3^ , T) (^c)) being D(c) a subbundle of D is maximally reducible 

to a driftless system {M,D,U) if 

(i) {M,D,U) is a kinematic reduction of {D , , 3^ , "D (^c)) , o-nd 

(ii) for every solution {'-y(t),u(t)) of Equations (5.2) satisfying '^{0) G 



{toO7)(0) 



, there exists a controlu G U such that {TD°"y{t),u{t)) 



is a solution of Equation (5.3). 



The characterization of maximally reducible mechanical control sys- 
tems defined on skew-symmetric algebroids is given by the following 
result. This a a generalization of the notion of being maximally re- 
ducible systems proved in |3] on Riemannian manifolds. 

Theorem 5.14. Let (D, S^, 5", 2)(c)) be a mechanical control system 
such that D(c) has locally constant rank and 3^ = 0. This mechanical 
control system is maximally reducible to a driftless system {M,D, U) if 
and only if the following two conditions hold: 



22 



BARBERO, DE LEON, MARRERO, MARTIN DE DIEGO, AND MUNOZ 



(ii) D(c) is geodesically invariant, that is, Sym*-°^''(D(c)) = I'(c)j 

where Sym^°°''(2)(c)) is the smallest distribution containing !D(c) and 
closed under the symmetric product (■ : ■)gi?. 



Proof. It follows the same lines as the proof of Theorem 8.27 in [1]. 

Suppose that the rank of 2) C D is k'. Note that X G ^{td) is a 
section of D if and only if the vertical lift of X restricted to D is tangent 
to 2). 

Assume that {D, S^, 0, 1)(c)) is maximally reducible to a driftless sys- 
tem (M,2),f7). 

If7:/C]R— j-Disa geodesic for the bundle metric 9^ with initial 
condition 7(0) G 2)o-(o), then 7 is a solution of Equations (5.2) with zero 



controls. Thus, by hypothesis, there exist controls : / — )■ f/ such that 
k' 

7(t) = ^u"{t)Xaia{t)), where a = td o ^ : I C R ^ M. Therefore, 
0=1 

7(1) C 2) and we can conclude that D is geodesically invariant. 
It remains to prove that D = 

Remember that D(c) = span{Yi, . . . , Y^}- Then it is sufficient to 
prove that Yg is a section of D for s = 1, . . . ,k in order to obtain 
D(^c) C. D. Having in mind that ^go is the geodesic spray on D and F/ 
is the vertical lift of Yg G TItd^^^) to D(c) , it is clear that the integral 



curves of ^gc + Y^ are solutions of Equations (5.2). By assumption, 
^gD + Y^ restricted to D is tangent to D. On the other hand, D is 
geodesically invariant, what implies that ^go restricted to D is tangent 
to D. Hence, (Y^)\^ is tangent to D and Yg is a section of D. 

Next, we will see that D C D(c). Let a be a vector in D. We consider 
the geodesic 7 : / C M — )• D with initial condition 7(0) = a. Hence, 
7(1) C 2) because 2) is geodesically invariant. Now, take the curve 
7 : / C M — )■ 2) given by 

7W=t7(^), VtG/. 
Note that r/j o 7 = r/j o 7 o r, where r : M — t- M is the map defined by 



r(t) = t^/2, and 7(0) = 0. From (3.1) it is straightforward to prove 
that 



f) (0) 



V?;)7 ) (0) = a. 
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On the other hand, define a curve a on M given by a = r^) o 7 such 
that 

k' k' 
7(t) = 5^n"(t)X„((r^o7)(t)) = ^"(t)X„((a o r)(t)). 

a=l 0=1 

Then, it is clear that 

k' 

J2^-it)pniXM^ o r)(t)) = pnilit)) = ° r) = |(r^ o 7), 



what imphes that {a o r, -u) is a solution of Equations (5.3). Therefore, 
there exist controls u such that (7,^) is a solution of Equations (5.2). 
In particular, 

« = V?"o)7(0)gD(c). 

The other implication is proved as follows: First, we prove that 
(M, 2), F) is a kinematic reduction of {D, S^, 3", 'D(c))- Let a : / C M ^ 
M be a curve such that it satisfies a{t) = Yl'l=i^°'i^)PD{Xa)io'{t)), 



i.e. Equation (5.3). Take 7(t) = Ylt=i^°'i^)-^aio'{t)) and note that 



7(1) CD and o 7 = a. Moreover, 

k' 

Thus 7 is pD-admissible. 



Now we have to prove that V^(.^-)7()f:) G 2)(c)(cr(t)). By hypothesis. 



2) = I'(c) and D(c) is geodesically invariant. Then by Theorem 3.2 



we 



conclude that V^/^^7(^) G D(c)(cr(t)) for all tel. Hence, there exist 



controls such that 7 is a solution of Equation (5.2). 



It remains to prove condition (ii) in Definition 5.13 Let (7,'u): / — ?■ 



D X U he a solution to Equation ( |5.3 ), then Q(c)(VLj)7(^)) = 0, that 



IS, 

V?",)7We%)Kt)) = 2)(a(t)) 

where cr = r/) o 7. 

By assumption 7(0) G !l)(cr(0)) and 7 is an integral curve of ^go +X^, 
where ^go is the geodesic spray associated with and X G r(rx)). 

Note that XJ^ is tangent to D and (^go)|^ is also tangent to D 
because D is geodesically invariant. Hence, if 7(0) G S)(cr(0)), then 
the integral curve of ^gu + X^ with initial condition 7(0) is entirely 
contained in D, that is, •yit) G D{(7(t)) for all t G /. 
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Then there exist u: I U such that 7(t) = Yl!L=i'^°'{^)^oi{<^{t)) and 

k' 



a=l 



Using the fact that 7 is p/j-admissible it immediately follows that a{t) = 
Eli«"(^)PD(X„)(a(t)) for all tel. □ 

Let us provide some results to find decoupling sections for the me- 
chanical control systems under consideration. 

Corollary 5.15. If a mechanical control system (D,S^ ,3^,'D(c)) with 
3^ = is maximally reducible, then all control sections are decoupling. 



Proof. This follows immediately from Theorem 5.14 and the definition 



of decoupling sections. □ 

The converse is not necessarily true. But it is true when D(c) has 
locally constant rank equal to one. 

Corollary 5.16. Let (D, S^, 5", !D(c)) be a mechanical control system 
where 2)(c) has rank one and 5" = 0. The following statements are 
equivalent: 

(i) There exist decoupling sections for {D, S^, 5", '^(c))- 
(a) The mechanical control system {D, Q^, 3^, '^(c)) is maximally re- 
ducible to a driftless system defined by decoupling sections. 



Proof. It is straightforward from Theorem 5.14, □ 



5.3. Examples. The first two examples are specific cases of Exam- 
ple |3.4[ As is shown, the example in Section |5.3.1 is not maximally re- 



ducible but admits rank-one kinematic reductions. Particular solutions 



to Hamilton-Jacobi differential equation are found using Theorem 5.9 



and Corollary 5.12 In Section 5.3.2 we compute particular solutions 



to Hamilton-Jacobi differential equation for a maximally reducible sys- 
tem. The snakeboard described in Section r5.3.3l has nonholonomic con- 
straints. Hence the use of skew-symmetric algebroids to find solutions 
to Hamilton-Jacobi differential equations is very natural. 

5.3.1. Planar rigid body with a variable- direction thruster. We refer 
to lU Section 7.4.2] for a detailed description of the system. The con- 
figuration space is M = X Consider local coordinates {x,y,6). 
Here the distribution D is the entire tangent space TM and = 0. We 



are in the case explained in Example 3.4 
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The riemannian metric is 

g™ = JdO ®de + m(da; ® dx + dy ® dy). 

The control vector fields in 1^(c) are 

cos6' d sinO d sin6' 9 cos^ d h d 

Y-^ = _l_ Y2 = -|- 

m dx m dy^ m dx m dy J 89 

In this example 

^i) = ^P""^- {^''- =-'''^lx^ ^ly^lle 

Note that {11,1^2,^3} is a 9"^^^-orthogonal basis. The skew-symmetric 
algebroid structure is defined as follows D = TM, |-,]rAf is the usual 
Lie bracket and ptm = Id^Af- In our adapted basis we have 

h 

{Yi, ^21tm = [^1, ^2] = r^^2 + 7^— r^^3, 

The non-zero Christoffel symbols for the associated Levi-Civita V^^^^ 
connection are 

i22- ^, i23- J , i32- ^, 

pi p2 h _ J 

^ 33 - 7, ) J- 21 - T , 31 



h' J + mh^' h{J + mh^)' 



h? ^„ h 

r3 p3 _ 



JiJ + mK^y J + mh^ 

Let us compute the symmetric products of the control vector fields 

<Yi:Yi>^n =0, <Yi:Y2>sn=Tl,Y2 + TlYs, 
<Y2:Y2><iD =2Tl2Yi. 

By Corollary |5.8[ Yi and Y2 are decoupling vector fields. However, 
!D(c) = span{Fi,F2} is not geodesically invariant because Sym''^-'D(c) ^ 
X'(c). Then according to Theorem 5.14 the mechanical control system is 
not maximally reducible to a driftless control system. As the involutive 
closure of the decoupling vector fields Yi and Y2 has maximum rank, the 
system is kinematic controllable, see Definition 5^ Then the motion 
planning is feasible. 
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Let US see if Yi and Y2 satisfy condition 
is, = g™(vf X,X) for all Y G Df. = spanlFg}- 



in Corollary 5.12 that 



s™(vr 



rLS™(r2,Fi) + rLs™(i3,ri) 



0, 



since {^1,^2 5^3} is a S"^^^-orthogonal basis. Thus Yi and I2 are both 
solutions to Hamilton- Jacobi differential equation because of Corol- 



lary 5.12 



Let us try to find more vector fields solution to Hamilton- Jacobi 
differential equation. In order to do this we have to find functions 



/ G C°°(M) such that /Yi and fY2 also satisfy condition (ii) in Corol- 



lary 5.12 According to Proposition 5.11 / must satisfy 



Y,if) = ^ --^-d- + cos9- + -- = ^. 

For instance, f{x,y,6) = g{x — hcos6,y — hsinO) where g: M."^ — > M 
satisfies the above partial differential equation. Hence, all vector fields 
fYi and /F2 are solutions to Hamilton- Jacobi differential equation, but 
not necessarily their linear combinations. 

So far we have found some particular solutions of Hamilton- Jacobi 
differential equation. Let us consider now the most general section in 
T{td), X = aiYi + a2Y2 + aala. Let us find functions ai,a2,a3 such 
that condition zx(d^(bgD(X)))(y) = for all Y G Vf^^ in Theorem 
is satisfied, that is, 

= - g^(ri, Yi){aiYs{ai) + a2air^2 + c^sc^i'^ls) 



5.9 



9''iY2,Y2)ia2Y,{a2) + a,a2T 



31) 



+ g^(r3, r3)(airi(a3) + a2>^2(«3) + «ia2r^i). 

It can be proved that this condition and Hamilton- Jacobi differential 



equation (5.4) are satisified, for instance, by 

g{x — h cos6, y — h sin^)y3, g{x — h cosO, y — h sin^^)y2 + ^3^3, 
g{x — hcos9,y — hsin9)Yi, /((x — hcos9)/h)Y2 + a^Y^. 

5.3.2. Robotic leg. We refer to [H Section 7.4.1] for a detailed descrip- 
tion of this system. The configuration manifold is M = IR+ x §^ x §^ 
with local coordinates {r,6,ip). The riemannian metric for the system 
is 

g™ = r7i(dr O dr + r^dO ® dO) + Jdifj dip, 



KINEMATIC REDUCTION AND THE HAMILTON- JACOBI EQUATION 27 

where m is the mass of the particle on the end of the extensible leg and 
J is the moment of inertia of the base rigid body about the pivot point. 

The control vector fields that span !D(c) are 

ri = — ^2 = --. 

mr^ d6 J dip ' m dr 

There are no constraints on the system, then as in the previous example 
D = TM, I , -Jtm is the usual Lie bracket and ptm = Wtm- 
The S ^-orthogonal distribution to D(c) is spanned by 

^ dijj de' 

Note that {Yi,y2;^3} is a S'^-orthogonal basis. From now on we con- 
sider coordinates adapted to this basis. Then the non-zero Christoffel 
symbols are: 

p2 ^ pi ^ J p3 ^ \ 

11 mr^' mr(J + mr'^)^ ■mr(J + mr^)^ 

r2 = = - , = ^ , 

13 21 mr{J + mr'^)^ '^'^ mr(J + mr'^)^ 

Fl _^IL_ r3 ^ r2 __1 

^23- j^^^2' ^23- ^^^^2' ^31- 

Fl - "IL F3 _ L f2 _ _ 

^32 -^^^^2' ^32 -^^^^2' ^33- ™- 

From here, it is easy to compute the following symmetric products: 
(Fi: Fi)gz, = ^F2, {Yi: ¥2)^0 = {Y2: ¥2)30 = 0. 



According to Theorem |5.14| the system is maximally reducible. Thus 
all the C°°(M)-linear combination of Yi and Y2 are decoupling sections 
and we can apply Corollary 5.12 to identify those decoupling sections 
that are solutions to Hamilton- Jacobi differential equation: 

= g^(VS^ (aiFi + «2>2), aiYi + a2Y2) 
= g^(yi, Y,) (aiFsK) + «ia2r^2) + S^(r2, ¥2) {a2Ys{a2) + aittsr^) 

It can be proved that if ai{r, 6, ip) = fi{r, 6 — ip) either for i = 1 01 i = 2 
and smooth functions /j: x S-*^ — )■ M, then the section aiYi -|- ^2^2 
is a solution to Hamilton- Jacobi differential equation. 



By Theorem |5.9[ we can also check that 

f{r,,p-e)Y, + a3Y,, f{r,,p-e)Y^ 
are solutions to Hamilton- Jacobi differential equation. 
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5.3.3. Snakeboard. We refer to Section 13.4] and [21] for a detailed 
description of this system. The configuration manifold is M = SE{2) x 
§^ X with local coordinates {x,y,9,ip,(p)- Consider the following 
physical parameters: the mass rric of coupler, mass rrir of rotor, mass 
m^t, of each wheel assembly, inertia Jc of coupler about center of mass, 
inertia Jr of rotor about center of mass, inertia of wheel assembly 
about center of mass, distance / from coupler center of mass to wheel 
assembly. 




The Riemannian metric in the given coordinates is 

g = {mc+mr+2m^){dxiS)dx+dyiS)dy) + {Jc+Jr+2{J^^+mJ^))d90d9 
+ Jrdip ^dij + 2J^d(j) (g) d0 + Jr{de (g) dtp + dip (g) dO). 

The constraints require that the wheels of the snakeboard roll without 
slipping. These constraints define a nonholonomic distribution D. As 
computed in Lemma 13.11] a S-orthogonal basis for D on Md = 
{{x,y,e,iP,^)\<P^±^} is 

f d d\ d d 

Xi = I cos (f) cos6— — h sin 6— — sin 0— = / cos (f)Vi — sin ( 



dx dy J 86 89^ 



where 



d_ 
50' 



, , , Jr-l COS (t) sin (t) , , Jr sin^ i 
a(0) = 



Ci(0) ' Ci(0) 

and 

Ci(0) = {rric + rrir + 2m^)Z^ cos^ + ( Jc + + 2( + m^l'^)) sin^ 0. 

Let us define a skew-symmetric algebroid structure on the vector bundle 
td'- D ^ Ml). If {61,62,63} is a local basis of T^td), then 

Pd(6i) = Xi, Pd(62) = X2, poies) = X3, 
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and 

[ei,ei]D = [62, 62]^= [e3,e3lD = 0, 
Iei,e2]D = 0, 

9([Xi,X3],Xi) , 9([Xi,X3],X2) 

= 9{X,,X,) S(X2,X2) 

S([X2,X3],Xi) , 9([X2,X3],X2) 

Hence the non- vanishing local structure functions for the bracket [■, -Jd 
are 

pi p2 pi p2 
^13^ *^13' '-'23' *^23- 

The bundle metric on the skew-symmetric algebroid is given by = 

We can construct the Levi- 



3.5 



S\dxmjjD as explained in Example 

Civita connection V^^ associated to the bundle metric having in 
mind its properties: 

^BC = 2^^^ {dCC^^AB + QAA^bc + aBB^Ac) 

for a S^-orthogonal basis, where g"^^ = I/qaa- 
The non- vanishing Christoffel symbols are 

Tis = oS' 5'llC23 + 2^13' ^23=2^' 5'22(0) + C23' 

p2 _ I„22^ pi lp2 r2 — lr,22^/ 

J-31-25' fl'll'-23 ~ 2'-13' ^32—29 922W)^ 

Tl,=gy^el,, ^ ^ Tl, = -l£'{g22el, + g^^el,), 

^21 = ~2^ (5'llC23 + ^'22613), r22 = 29 922^^2- 

The mechanical control system on the skew-symmetric algebroid is given 
by (D,g^,0,X)(e)), where 

^ _ [ 1 



where C2(0) = (mc -|- m.^ -|- 2171^)1^ cos^ -|- ( Jc + 2(7^^, -|- rrii^l^)) sin^ 

The control forces are torques, one actuating the rotor and the other 
one actuating the wheels. It is easy to prove that 
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In order to find solutions to Hamilton- Jacobi differential equation, let 
us first check if there exist any decoupling section of r(!D(c)) so that we 



can use Corollary 5.12 



(X2: X2)c,D = 2TI2X3 e D(c), 

(X3:X3)gD=0eD(,). 

Hence both control sections are decoupling. As in the previous example 
the mechanical control system is not maximally reducible to a driftless 
system because 

{X2 : X3)gD = (r23 + r32)Xi + (123 + ^32)^2 4- "^{c)- 

Let us check if the decoupling sections satisfy condition (ii)| in Corol- 



lary 5.12 



r?2g^(X3,X2) = o, 
g^(vi;x3,X3) = ri3g^(Xi,X3) + r23g^(X2,X3) = o 



because {Xi,X2,X3} is S -^-orthogonal basis. Thus X2 and X3 are both 
solutions to Hamilton- Jacobi differential equation because of Corol- 
lary 5.12 Let us see if these solutions can generate more decoupling 
sections being solution to Hamilton- Jacobi differential equation. As 
proved in Proposition 5.11, /X2 and /X3 with / G C°°(M) are solu- 
tions to Hamilton- Jacobi differential equation if and only if 



cos U— — h sm U— 
ox oy 



sm( 



de 



0. 



2xsm.(i) ^ , . n ,cos0 _ , , 
+ 2cos0sm6',?/ — L— — - cos 6, ip,(p 



I 



sm I 



Xi(/) = ^ /cos 
If we take any function 

f{x,y,e,ip,(f)) = g 

where 5?: — )■ M is a solution to the partial differential equation 
Xi{f) = 0, then /X2 and /X3 are also solutions to Hamilton- Jacobi 
differential equation. 

Consider now a general section aiXi + 0:2X2 + 0:3X3 i n T{t £)). Con- 
becomes 



dition ^xd^(bgD(X))(r) = for all Y G Vf^^ in Theorem 



5.9 



g^(Xi,Xi)(02X2(0i) + 03X3(01) + 030ir3i + 0203123 + 0302132) 



g^(X2,X2)(02Xi(02) +03021^3) 

g'^(X3,X3)(o3Xi(o3) + aiasTl-^ + 03021 



12 J 



0. 



For instance, for any Oi G C°^(M) OiXi satisfies this condition. It will 
be a solution to Hamilton- Jacobi differential equation if it fulfills (5.4), 
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that is, 

=g^(VS"^^aiXi,Xi) = aiXi(ai)g^(Xi,Xi) + alVl^S'^iXs, X,) 
=«iXi(ai)g^(Xi,Xi). 

Equivalently, Xi(ai) = 0. Thus aiXi is also a solution to Hamilton- 
Jacobi differential equation for any 

/'2xsin0 COS0 

ai[x, y, 0,%j,(p) = g { ; h 2 cos^smb', y — I cos 0, ip, 

\ I smcp 

where 51 : M. 



6. Future work 



The results in this paper extend the notion of decouphng sections 
for mechanical systems with nonzero potential. A decoupling section 
for those mechanical systems is a section that satisfies the assumption 
and condition (ii) in Theorem 5.9 The future research line consists of 
taking advatange of this geometric description of decoupling sections 
to do motion planning for mechanical systems with nonzero potential. 
One of the key points to succeed in motion planning is that not any 



reparametrization of decoupling sections in the sense of Theorem 5.9 
is again decoupling in the sense of Theorem |5.9[ In order to define 



the suitable notion of reparametrization Proposition 5.11 seems to be 
useful. 
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